Orbital degeneracy removed by charge order in triangular antiferromagnet AgNiC>2 
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We report a high- resolution neutron diffraction study on the orbitally-degenerate spin- 1/2 hexag- 
onal metallic antiferromagnet AgNiCb- A structural transition to a tripled unit cell with expanded 
and contracted Ni06 octahedra indicates y/3 x y/3 charge order on the Ni triangular lattice. This 
suggests charge order as a possible mechanism of lifting the orbital degeneracy in the presence 
of charge fluctuations, as an alternative to the more usual Jahn- Teller distortions. A novel mag- 
netic ground state is observed at low temperatures with the electron-rich S = 1 Ni sites arranged 
in alternating ferromagnetic rows on a triangular lattice, surrounded by a honeycomb network of 
non-magnetic and metallic Ni ions. We also report first-principles band-structure calculations that 
explain microscopically the origin of these phenomena. 

PACS numbers: 75.25,+z, 71.45.Lr 



Electronic systems on triangular lattices have at- 
tracted wide interest theoretically and experimentally 
due to the possible existence of unconventional ground 
states stabilized by the frustrated lattice geometry [l|. 
The layered cobaltate NaCo02 with a triangular lattice 
of orbit ally- nondegenerate Co 3+ ions [2] shows a num- 
ber of unusual phases upon hole doping via deintercala- 
tion of Na, as well as superconductivity upon hydration 
at a particular doping [1]. The addition of orbital de- 
generacy creates more instabilities for cooperative order 
and the behavior in the presence of metallic conductivity 
is not well understood and largely experimentally unex- 
plored. In more localized systems such as NaNi02 with 
spin- 1/2 Ni 3+ ions in the t% g e\ configuration the two- 
fold e g orbital degeneracy is lifted by a Jahn- Teller (JT) 
structural distortion to a monoclinic structure [3] , which 
splits the e g band and opens a gap. However it is un- 
clear if such a mechanism applies in less localized sys- 
tems where the local tendency for JT distortions com- 
petes with the charge transfer between sites, allowed by 
an efficient metallic screening [4|. The metallic layered 
silver nickelates AgNiC>2 [H, @] and Ag2Ni02 [7] are ideal 
candidates to explore this physics. 

Here we report high-resolution structural studies of 
AgNi02 and find that no JT distortions occur but the 
ideal structure is distorted to a hexagonal tripled unit 
cell with expanded and contracted Ni06 octahedra. This 
implies that the orbital degeneracy of the formally Ni 3+ 
sites is lifted by a novel \/3 x \/3 charge ordering (CO) 
that occurs via the charge transfer 3e g — > e 2 g + e g - 5 + e^ 5 , 
which to the best of our knowledge is the first exper- 
imental example of such a CO pattern in a triangular 
system. Using band-structure calculations we show that 



this CO is a 2D analogue of the metal-insulator transi- 
tion in the 3D KH'iO^ [8[ attributed to the charge transfer 
2e g — > e 2 g + e® [4]. However, in contrast to the 3D nick- 
elates in AgNi02 the triangular lattice geometry allows 
the system to remain metallic in the CO phase, as itiner- 
ant electrons can hop on the inter-connected honeycombe 
network formed by the two electron-depleted sublattices 
(e^ 5 ); those sites remain magnetically unordered at base 
temperature. The remaining one third of Ni sites are 
strongly magnetic (e g ) and form an effectively tripled an- 
tiferromagnetic (AFM) triangular lattice, which orders in 
an unusual collinear 2x1 structure of alternating stripes. 

Powder samples of 2H-AgNi02 (a polymorph [5[ of 
the better known rhombohedral 3R-polytype [6]) were 
prepared as described in Ref. [5] and neutron diffrac- 
tion indicates a nearly pure hexagonal phase with less 
than 1% admixture of the rhombohedral poly type. Neu- 
tron diffraction patterns at 300 K were collected using 
the high-resolution time-of-flight diffractometers OSIRIS 
and HRPD at the ISIS pulsed neutron source and rep- 
resentative results are shown in Fig. yielding the 
hexagonal space group P63 /mmc as proposed earlier 
[5[, with triangular Ni, O, and Ag layers stacked as 
ABBBA (Ni-O-Ag-O-Ni). However the neutron data 
shows additional diffraction peaks [see Fig.Q}))-d)] which 
can be consistently indexed by an ordering wavevector 
&o = (1/3,1/3,0) in the original undistorted structure, 
indicating a tripled unit cell in the hexagonal plane. To 
test the accuracy of the indexing we have refined the lat- 
tice parameters using only the main peaks and only the 
supercell peaks (more than 20 observed) and obtained 
the same values to within the experimental accuracy of 
10 _4 A. Furthermore, upon cooling to 2 K the super- 
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cell peaks were displaced in Q following the lattice con- 
traction but maintained their commensurate index with 
respect to the main peaks. This further confirms that 
the supercell peaks are due to a modulation of the main 
structure rather than an extra impurity phase. 




FIG. 1: a) (color online) 300 K neutron powder diffrac- 
tion pattern and b)-d) zoomed-in regions showing some of 
the prominent supercell peaks which indicate a distortion of 
the high-symmetry P63 /mmc structure where all Ni sites are 
equivalent. e)-g) Supercell peaks are not observed in xray 
data indicating that the distortion mainly involves displace- 
ments of oxygen ions. Solid lines in all panels are fits to the 
model with periodic arrangement of expanded and contracted 
Ni06 octahedra shown in Fig. [2] and peak indexing is in the 
triple unit cell (space group P6322). 



X-ray data collected using a Phillips Cu K a powder 
diffractometer [Fig Hfe)-g)] do not show supercell peaks, 
indicating that the distortion involves mainly the light 
oxygen ions. Thus we refined the neutron data in a 
y/3 x y/3 hexagonal unit cell assuming only oxygen dis- 
placements. We eliminated the space groups that were 
not compatible with the observed low-temperature mag- 
netic structure (to be described later) dictating at least 
two different Ni crystallographic sites. The highest sym- 
metry space group in which both the structural and mag- 
netic data could be described is P6322 (#182). Good 
agreement with the data is obtained (solid line fits [9[ in 
Fig [1]) and the extracted unit cell parameters are listed 
in Table H As a further check of the space group identi- 
fication we have performed full structural optimization 
calculations (as described later) starting from a lower 



TABLE I: Unit cell lattice parameters in the ideal (P63/rarac) 
and CO (P6322) structures at 300 K. Internal parameters 8, 
C, f, and e are 0, 0.08050(5), and 0.0133(2) (exp.) and 0, 
0.0792, and 0.0102 (calc), respectively. 



P6 3 /mmc (no. 194) 



a = 2.93919(5) A 
c = 12.2498(1) A 



P6 3 22 (no. 182) 



a = 5.0908(1) A 
c = 12.2498(1) A 
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symmetry structure (P63 space group) and those invari- 
ably converged to the higher P6322 symmetry, with the 
optimized positions agreeing well with the experimental 
ones (Table U) [To|. 
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FIG. 2: (color online) Bottom MO2 layer showing the ex- 
panded Nil06 octahedra (large circle) surrounded by a hon- 
eycomb network (thick lines) of contracted Ni2 and Ni3 sites 
(small circles). The small balls are oxygen ions and the small 
arrows indicate the displacements from the ideal structure. 

The resulting structure in one MO2 layer is shown in 
Fig. [2j There are three inequivalent Ni sites. Oxygens 
breathe in towards Ni2 and Ni3 ions (c?M2-o = <^Ni3-o = 
1.934 A) and out from Nil (d N ii-o = 2.022 A). No Jahn- 
Teller distortion is present. The different Ni-0 distances 
indicate that Nil is electron rich (Ni 2+ in a CO scenario) 
and Ni2,3 are electron depleted (Ni 3,5+ ). In the adjacent 
MO2 layer Nil and Ni3 sites swap places, leading to a 
zigzag arrangement of the (expanded) electron-rich Nil 
sites along c [111 ]. 

Earlier susceptibility data [5] showed dominant an- 
tiferromagnetic couplings with a Curie- Weiss constant 
= -107.6 K and AFM order below Tat=22 K. Below 
this temperature additional peaks appear in the neutron 
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FIG. 3: a) (color online) Magnetic diffraction pattern (dif- 
ference 4 K - 300 K) indexed by wavevector k = (1/2,0,0). 
Solid line is a fit to the spin order depicted in b): In one layer 
ordered spins (Nil) form alternating ferromagnetic rows, zb 
symbols indicate spin projection along the c-axis and gray 
balls are unordered sites (Ni2 and Ni3). Diamond-shaped 
box is the unit cell, c) 3D view of the magnetic order along 
c-axis. Band-structure calculations indicate that the effec- 
tive FM inter-layer coupling between the strongly-magnetic 
Nil spins (large arrows) occurs via AFM couplings to a very 
small moment at the Ni3 site (small white arrow). Inset in 
a) shows the 2D hexagonal Brillouin zone with locations of 
magnetic Bragg peaks (filled stars for the stripe domain in 
(b), open stars for equivalent domains rotated by ±60°). 



diffraction pattern, see Fig.[3fci). These peaks are indexed 
by the ordering wavevector k — (1/2,0,0), as shown in 
Fig. [3^) inset. The best fit was obtained when only one Ni 
sublattice was ordered, either Nil or Ni3, with a moment 
of 1.522(7) hb along the c-axis. As discussed above, this 
must be the strongly-magnetic Nil (Ni 2+ ) with a large 
available moment of 5=1 (2/i#), because Ni3 and Ni2 are 
Ni 3 5+ with only a small available moment. 

The observed magnetic structure is shown in Fig. 03d): 
in each layer ordered spins are arranged in alternating 
FM stripes, in the adjacent layer stripes are parallel but 
have an in-plane offset that follows the c-axis zig-zag ar- 
rangement of the electron-rich Nil sites, see Fig. [3h). 
This magnetic structure implies that (a) the exchange 
between nearest in-plane magnetic neighbors is AFM, (b) 
there is an easy-axis anisotropy —DS^ and (c) the net 
interplanar interaction with the 3 magnetic neighbors in 
the layer above (below) [see Fig. [3b)] is FM. 

To gain more insight into the physical origin for the 
observed charge and magnetic orderings, we have per- 
formed density functional (DFT) calculations 1_2|, 131 . 
Since Ag x Ni02 is metallic and only weakly correlated 14j, 



introducing a Hubbard U in the calculation is not nec- 
essary. Optimizing the oxygen positions in the y/3 x y/3 
unit cell precipitates the CO and reproduces the observed 
structure in Fig. [21 indicating that the driving force for 
charge disproportionation is entirely accounted for in the 
LDA. The calculated total magnetization for the three 
Ni ions varies between 1.3 and 1.5 depending on the 
exact O position, and is nearly entirely located at the 
expanded Nil site in good quantitative accord with the 
neutron data. The computed ferromagnetic band struc- 
ture (Fig.SJ) shows fully exchange-split Nil derived bands 
with the magnetic moment diminished from its formal 
count of 2 hb by hybridization. One can also see that 
Ni2 and Ni3 states are not visibly exchange split. The 
Ag sp band is entirely above the Fermi level. This elec- 
tronic structure can be visualized as a magnetic insula- 
tor formed by Ni 2+ (cf. NiO) with a strong tendency 
to magnetic order, superimposed on a Ni 3 5 metal. The 
latter has larger bandwidths, due to the smaller Ni-0 
distance, and by itself does not satisfy the Stoner crite- 
rion for metallic magnetism, IN(Ep) > 1, where / is the 
atomic Stoner parameter. 

The natural tendency for charge order in a MO2 layer 
can be explained by energetic arguments. The Hund's 
rule energy gain in a metal in DFT is M 2 (I — 7V _1 )/4, 
where N is the average density of states per spin. With- 
out charge ordering, Mi = M 2 = M3 = ljiB and 
Ni = N 2 = N 3 = N. After charge ordering, Mi « 2/i B , 
M 2 = M 3 « and JVi > N > N 2 = N 3 . The CO state is 
energetically favored since 4(/-7V 1 _1 )/4 > 3(/-7V _1 )/4. 
In insulators such as NaNi02 the Hubbard repulsion 
would forbid this mechanism, but in this metallic sys- 
tem it is well screened and cannot prevent CO [15]. Our 
band-structure calculations reproduce fully the observed 
3D CO pattern, both the in-plane \/3 x y/3 order as well 
as the zig-zag arrangement of the electron rich sites along 
c. We also found a metastable solution where the CO lay- 
ers form straight columns along c, i.e. the electron rich 
and strongly magnetic Ni 2+ occupy the Ni2 site; this 
indicates that the in-plane CO is very robust and inde- 
pendent of the inter-layer magnetic interactions. 

To study the stability of the observed magnetic struc- 
ture we have calculated ground state energies for a num- 
ber of potential structures using DFT for the full mag- 
netic unit cell (12 formula units and lower symmetry, C2). 
In addition to the observed stripe-order with the zig-zag 
ferromagnetic pattern along c shown in Figs.[3)3-c) (called 
stripe-F), we also found stable solutions for hypothetical 
ferromagnetic (FM), A- type AFM (alternating FM lay- 
ers), as well as in-plane stripe-order but with a zig-zag 
AFM pattern along c (stripe-A, every other plane flipped 
compared to stripe-F). The results are listed in Table OH 
the state found experimentally has the lowest energy. 

The magnetic structure indicates that the in-plane ex- 
change between magnetic neighbors is AFM. This arises 
naturally because Nil 2+ sites are insulating and lack 
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TABLE II: Energies of different magnetic configurations for 
2H-AgNi02- All energies are expressed per magnetic Nil ion 
and as a difference from the calculated (and observed) mag- 
netic ground state. 





stripe-F stripe-A FM A-type AFM 


Energy/Nil (meV) 


7 23 21 



both metallic Stoner ferromagnetism and 90° FM su- 
perexchange, but do have a classical AFM superexchange 
via the Nil-O-O-Nil path, with effective hopping t e ff ~ 
tp da t pp7r /(E ( i — E p ) 2 . The appearance of collinear order 
instead of the conventional 120° order for a triangular 
AFM implies a significant easy-axis anisotropy. This is 
normally provided by spin-orbit induced coupling to the 
crystal field, but the Nil ion is close to a t^e 2 configu- 
ration with zero orbital moment. Further experimental 
and theoretical studies of the single-site anisotropy in 
this system are under way. The fact that the effective 
interplanar Nil-Nil exchange is FM is rather surprising 
since there is a priori no apparent physical mechanism 
for FM couplings. Indeed, the magnetic species, Nil, 
forms insulating bands and Ag is nonmagnetic and ab- 
sent from the Fermi energy, eliminating double exchange 
as a possibility. Ni2 and Ni3 form quasi-2D bands with 
small potential for sizeable RKKY interaction, and mag- 
netic impurities in the Ag layer that could mediate a FM 
exchange have been excluded experimentally. 




FIG. 4: (color online) Ferromagnetic band structure of 2H- 
AgNi02 along symmetry directions in the Brillouin zone. 
Blue/red (light/dark shades) correspond to the up/down spin 
directions and circles indicate the relative weight of the Nil 
states. Inset: Partial densities of states for the three Ni 
sites in states/eV/triple cell. Top/bottom: spin majority 
(up) / minority (down) . 

Performing calculations for a wider range of parame- 
ters, we find that depending on the exact positions of 
Ag and O the interplanar interaction can change ampli- 
tude and even sign, and that an effective FM interlayer 
coupling can appear if the weakly-magnetic site Ni3 has 
a small ordered moment (0.10-0.15 Hb)- The strongest 
inter-planar interaction is AFM superexchange between 
the Ni3 and Nil directly above (it occurs via a near 
collinear O-Ag-0 path), whereas the inter-planar Nil-Nil 



exchange is much weaker. If the Ni3 moment is aligned 
antiparallel to the Nil above it then the structure in Fig. 
[3t) has a net energy gain as it also satisfies two out of the 
three in-plane AFM Ni3-Nil couplings. Our diffraction 
data in Fig. [3fci) cannot prove or exclude the existence of 
very small moments on Ni3, but the data is consistent 
with a Ni3 moment of ~ 0.1/i# anti-aligned with the Nil 
directly above it, as it comes out of the calculations. 

To summarize, we have reported high-resolution neu- 
tron diffraction in the orbitally-degenerate 2H-AgNi02 
which observe a periodic arrangement of expanded and 
contracted NiOe octahedra, naturally explained by a 
three-sublattice charge order pattern on the triangular 
lattice of Ni sites. We have proposed that due to Hund's 
rule coupling and metallic screening of the Hubbard U 
repulsion (and static screening due to oxygen displace- 
ments) the CO is the favored mechanism to lift the orbital 
degeneracy as opposed to the conventional Jahn- Teller 
distortions that occurs in more insulating systems. An 
unusual magnetic order is observed at low temperatures 
with only one third of sites carrying a magnetic moment 
with an unexpected collinear stripe order pattern on an 
antiferromagnetic triangular lattice. Our results indicate 
the complex cooperative charge and magnetic order pat- 
terns that can occur in orbitally-degenerate metallic sys- 
tems on frustrated lattices. 
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